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Abstract 

This work is devoted to the formulation and derivation of the a—K—fj,/ gamma distribution which 
corresponds to a physical fading model. The proposed distribution is composite and is constituted 
by the ct— k —p non-linear generalized multipath model and the gamma shadowing model. It also 
constitute the basis for deriving the a—n —p Extreme/ gamma model which accounts for non-linear 
severe multipath and shadowing effects and also includes the more widely known a —p and k—h 
models which includes as special cases the Rice, Weibull, Nakagami-m and Rayleigh distributions. The 
derived models provide accurate characterisation of the simultaneous occurrence of multipath fading and 
shadowing effects. This is achieved thanks to the remarkable flexibility of their named parameters which 
have been shown to render them capable of providing good fittings to experimental data associated with 
realistic communication scenarios. This is also evident by the fact that they include as special cases the 
widely known composite fading models such as the recently reported K—fj ,/gamma model and the novel 
a— /i/gamma model. Novel analytic expressions are derived for the corresponding probability density 
function of these distributions which are expressed in a convenient algebraic form and can be efficiently 
utilized in the derivation of numerous vital measures in investigations related to the analytic performance 
evaluation of digital communications over composite multipath/shadowing fading channels. 
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Index Terms 

a—n Distribution, k—/j, Distribution, k—// Extreme /gamma distribution, non-linear fading, shad¬ 
owing, composite fading channels, probability, performance analysis. 

I. Introduction 

It is widely known that fading is an effect that significantly degrades communication signals 
during wireless propagation. A common approach for accounting for this phenomenon has been 
possible through exploitation of suitable statistical distributions. To this effect, statistical models 
such as Rayleigh, Nagakami-m, Weibull and Nakagami-q (Hoyt) have been shown to be capable 
of modelling small-scale fading in non-line-of-sight (NLOS) communication scenarios whereas 
Nakagami-n (Rice) distribution has been typically utilized in characterizing multipath fading in 
line-of-sight (LOS) communication scenarios, [lU-EO] and the references therein. Capitalizing on 
these models, M. D. Yacoub proposed three generalised fading models, namely, the a— /x, the 
k~h, the r)~n models and subsequently the A—// and the K—fi Extreme models, lf4’l- [fT0l . These 
models along with the recently proposed a—n—fi Extreme fading model in fTO are useful thanks 
to the remarkable flexibility offered by their named parameters which render them capable of 
providing adequate fittings to results obtained by field measurements. Their usefulness is also 
evident by the fact that they include as special cases all the above small-scale fading distributions. 

It is recalled here that a fundamental principle of wireless radio propagation is that multipath 
and shadowing effects occur simultaneously. As a result, in spite of the undoubted usefulness of 
the aforementioned fading models, they all ultimately fail to account concurrently for both shad¬ 
owing and multipath fading. In other words, the characterisation offered by the aforementioned 
fading models is limited to modelling either the one or the other effect. Based on this crucial 
limitation, the need for composite statistical models that can provide adequate characterization 
of the fading effect as a whole became necessary |[ITl-|f3ll. 

Motivated by this, the authors in |fT2| proposed the Rayleigh/gamma fading model, which is 
also known as A'-distibution, (AT). In the same context, Shankar in lfT3l exploited the flexibility of 
Nakagami-m distribution, which includes Rayleigh distribution as a special case and introduced 
the Nakagami-m/gamma composite distribution - or generalised AT-distribution, {Kg). Likewise, 
the Weibull/gamma composite distribution was reported in lfl4ll while the k—/ i/gamma, the K—ji 
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Extreme/ gamma, the //—///gamma and the A—///gamma distributions were proposed in Ifl5l - 

m. 

The aim of this work is the derivation of novel analytic results for the a— k —///gamma 
composite fading model. Subsequently, by using this model as a basis, the a —///gamma and 
the a-n-fi Extreme /gamma composite fading distributions are additionally derived. Impor¬ 
tantly, unlike all existing composite fading models, this models are capable of accounting for 
the non-linearity of the wireless fading channel which is characterized by the parameter a. 
After formulating these models, novel analytic expressions are derived for their corresponding 
probability density function (pdf). The validity of the offered expressions is justified through 
comparisons with numerical results while their behaviour is assessed under different parametric 
scenarios. Importantly, owing to the relatively convenient algebraic representation of the offered 
expressions, they can be considered useful mathematical tools that can be efficiently utilized in 
studies related to the analytical performance evaluation of digital communications over non-linear 
multipath/shadowing and severe composite fading models. Indicatively, they can be exclusively 
employed in the derivation of algebraically complex analytic expressions for various performance 
measures such as error probability, channel capacity and higher order statistics, to name a few 

m-m. 

The remainder of this paper is organised as follows: Section II revisits the foundations of 
a—K—fi and the gamma distributions. Subsequently, Sections III, IV and V are devoted to the 
formulation and analytical derivation of the a— k— ///gamma, the a —///gamma and the q —//—// 
Extreme /gamma fading models, respectively. Finally, discussions on the potential applicability 
of the derived expressions in wireless communications along with closing remarks are given in 
Section VI. 


II. The a-K~n and Gamma Fading Distributions 
A. The a—K—fi Fading Model 

The ol — k — // distribution is a recently proposed fading model that accounts for the charac¬ 
terization of small-scale variations of mobile radio signals in LOS communication scenarios. 
It is written in terms of three physical parameters, namely a, k and //; the first one denotes 
the non-linearity parameter, // is related to the multipath clustering while n represents the ratio 
between the total power of the dominant components and the total power of the scattered waves. 
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For a fading signal with envelope R and r = “ \JE(R"), the a—n—p envelope pdf is given by 

m. 


r 
r 

where E(.) denotes expectation, f is the root-mean-square ( rms ) value of R and I u (x) represents 
the modified Bessel function of the first kind with argument x and order u If30l . For the case of 
normalised envelope P — R/r, equation ([]]) can be equivalently expressed as 


a(l+M) -i 

) 2 (l) 


Pr{t) = 


f( l + K)-^«^ e K"+^+«3S) 


Pp(p) = 
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while the corresponding cdf is expressed as, 


( 2 ) 


f>(p)=E 


/c*/iT (i + /i, /i(l + k)p“) 
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i\eP K r(i + /r) 

= 1 - (y/2/m, p% \/2p(l + k)J 


(3) 

(4) 


where T(a, x) and Q m (a, b) are the upper incomplete Gamma function and the generalized 
Marcum Q-function which are defined as [f30l . OT1 . Importantly, the a—n—ji fading model 
includes as specific cases the more widely known a—p and n—p generalized models ll27ll . In 
more details, for the specific case that a = 2, the a—n—p distribution reduces to the more widely 
known k— p fading model. Likewise, for the specific case a = 2 and k — 0, the Nakagami-m 
distribution is deduced while for m — 1 one obtains the Rayleigh distribution. In the same 
context, the Rice distribution can be obtained by setting n = k (where k denotes the Rice 
parameter) as well as for the specific case a = 2 and p — 1. Finally, the Weibull distribution 
can be obtained by setting k = 0 and p — 1. 

By setting pw{w) = pp{y/w)/2^/w according to [j8l, the corresponding power (pdf) is deduced 
in a straightforward manner, namely, 


p P (w) = 


apw ( 4 M) l Ifj- 1 ^2 p\J k(1 + n)w a /‘ 2 ^j 


2«V(1 + K yP^ e ^+™ a / 2 +™ a/2 ) 

In the same context, the corresponding statistical moments are given explicitly by, 


( 5 ) 
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e [/>'] = rM 


( 6 ) 


e^ K (l + K) l / a /j l / a 

where the notation i F\ (a: b; x) denotes the Kummer confluent hypergeometric function, 

Finally, it is recalled that the Nakagami parameter m is regarded as the inverse of the variance 
of the normalized power of the fading signal, i.e. m = Var ~ 1 (P 2 ), [(8]|. This term is linked with 
the k and // parameters through the following relationship, 


m = 


+( 1 + n) 2 

1 + 2 K 


(7) 


As already mentioned, the a—K—p fading model includes as a special case the a-n-ft 
Extreme distribution which subsequently includes the n—p Extreme model for a = 2. It is 
recalled that the validity of this model is based on that the occurrence of large number of 
paths is not typically the case in wireless radio propagation over enclosed environments such 
as aeroplanes, trains and buses ll28ll . [|29ll . The reason underlying this principle is that such 
environments are known to be characterised by only a few number of paths. As a consequence, 
severe fading conditions, even worse than Rayleigh, are ultimately constituted lUOll . ll28Tl . [[291 . 
The envelope pdf of the k —/i fading model was derived in llTOl and is given by, 


p P (p) = 4 me- 2m ( 1+p2 )/ 1 {Amp) + e~ 2m 5{p) (8) 

while the a—K—fi model which additionally includes the non-linearity of the wireless channel 
was recently proposed in fl8l and is defined as follows, 


For a 


Pp{p ) 


2amli (4mp Q//2 ) 

pl-f e 2m(l+p“) 


+ e~ 2m 5{p) 


2 equation © reduces to ®. 


(9) 


B. The Gamma Fading Model 

It is recalled that log-normal distribution has been regarded the optimum statistical model 
for characterising the shadowing effect, 0]|-[|3l. Nevertheless, in spite of its usefulness, it 
has been largely shown that when it is involved in combination with other elementary and/or 
special function, its algebraic representation becomes intractable. This is particularly the case 
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in studies related to the analytical derivation of critical performance measures in the field of 
digital communications over fading channels. Motivated by this, the authors in lfl2l proposed 
the gamma distribution as an accurate substitute to log-normal distribution. Mathematically, the 
envelope pdf of gamma distribution is given by lfT2l eq. (4)], namely, 


Py(v) 


y b 1 e n 
T(b)Q b ’ 


y> o 


( 10 ) 


where the term b > 0 is its shaping parameter and Q = E(Y 2 ). The gamma fading model 
has been shown to provide adequate fitting to experimental data that correspond to realistic 
fading conditions. In addition, it is evident that its algebraic representation is quite tractable and 
therefore, easy to handle both analytically and numerically. As a result, it has been undoubtedly 
useful in characterising shadowing and for this reason it has been exploited in the formulation 
of the K and K G composite multipath/shadowing models, [fT2|. Ifl3ll . 


III. The a- hi-fi /gamma Fading Distribution 
A. Model Formulation 

According to the basic principles of statistics, the envelope pdf of a composite statistical 
distribution is deduced by the superposition of two or more statistical distributions. In the 
present case, this is realized by superimposing one multipath and one shadowing distribution. 
Mathematically, this principle is expressed as, 


p x ( x )= p X \y( x I y)p Y (y) d y (U) 

Jo 

where p xlY (x \ y ) denotes the corresponding multipath distribution with mode y. Based on this 
principle, the corresponding n — ///gamma composite fading distribution is formulated by firstly 
setting in (Q]), r = x and r = y and then substituting in (flTT) along with equation (flOl) . To this 
end, it follows that, 


Px( x ) =A 




i i o(l+M) 

y 2 


-ftgMi+K) 


y_ 

en 


dy 


( 12 ) 


where 
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A = 


1 , \ 1 ±M 

CK/itC 2 X (1 + K) 2 

e^nfir T(b)tt b 


(13) 


By setting in (fl2l) . u = y a and thus, du/dy = ay a 1 and y = u 1 /", the above expression can 
be equivalently re-written as follows: 


Px( x ) 


A 1 (2/i-\/«(l + «) ^ 1 

— / --- ~A~du 

rv I -i 6 . M+l i a; a u*-/ a 


(14) 


B. A Novel Expression for the Envelope pdf 

It is evident from (fl4l) that the derivation of an analytic expression for the envelope pdf of 
the oi — k — p ,/gamma distribution is subject to evaluation of the involved integral. To this end, it 
is recalled the the I,fx) function can be equivalently represented in terms of the infinite series 
in [30, eq. (8.445)], and the polynomial approximation in 132'. eq. (19)], namely, 
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Therefore, by making the necessary change of variables 1 it immediately follows tha 
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+/ 
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Subsequently, by recalling that T(x) = (x — 1)! and substituting in ([T71) in (fT4l) yields, 


(17) 


Px\ x ) 


r (n + i)x a ^+ l ) 1 p^ + 2 l A(i + nY +l 
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n du 


(18) 


■"V"" 

Xi 


Importantly, the Xi integral in (fl8l) can be solved in closed-form. To this end, by expressing the 
corresponding exponentials in terms of the Meijer G-function according to |f33l eq. (8.4.3.1)] 
and |3_3, eq. (8.4.3.2)], yields, 


'The polynomial approximation is used in the present analysis. 

2 As n —1 oo, the polynomial approximation reduces to the infinite series. 
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By setting in (fl9l) . z = A" and thus, u = z a and dz/du = u<* 1 / a, one obtains, 


r<i,o 

-oo ^0,1 \ Q 


Xi — a 


0,1 


Ot i 
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J 0 z a(n+l)— 6+1 

In its current form, the above expression can be expressed in closed-form according to 
(2.24.1.1)] yielding, 


Xi = 


b-a(jj,+l )+1 

^0,1+Q! 

~L j 


n a a c 


1 -I 1 — b+a(p,+l) 
i ’ i > 7 


(19) 
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eq. 
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Therefore, by substituting (fUT) in (fl8l) . an analytic expression for the envelope pdf of the 
a—K—fi /gamma composite fading model is deduced, namely, 


p x ( x ) = 

1=0 


T (n + l)x a(ji+l ')~ 1 fj, fi+2l K l (l +k)^ 1 a b ~ a ^ +l)+1 o,i+a 
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’ Z 
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( 22 ) 


which to the best of the authors’ knowledge, it has not been previously reported in the open 
technical literature. 


C. A/7 analytic expression for the pdf of k— p/gamma model 

As already mentioned, analytic expressions for the pdf of the k— /i/gamma distribution was 
reported recently in iT5l . llT8ll . By recalling that the K—p distribution is a special case of the 
a—n—p distribution, an additional analytic expression for the pdf of k—/// gamma distribution 
can be obtained by setting a = 2 in (l22l) . Therefore, one obtains straightforwardly, 


A r(n + l)x 2 ^ +v, ~ l p^ 21 k 1 {1 + nY l+l 2 fc - 2 ^+9+5 0 3 / An 2 

Z!T(n — l + l)T(/i + Z)T(&)n 2Z_ A K y / 7rfZ 2 (i i +0 3,0 1^(1 + k)x 2 


1~Z>-|-2(/a-]-Z) I 

(23) 


The above expression could be useful in cases that the analytic expressions in [fl5ll . llT8l are 
analytically intractable since although the Meijer G-functions are considered rather laborious, 
they commonly allow the derivation of corresponding solutions in closed-form. 


























IV. The a - ji / gamma Composite Fading Model. 


A. Model Formulation 

Unlike reported results in generalised linear and non-linear multipath/shadowing fading mod¬ 
els, no results have been reported on composite models that are based on the a—p, fading 
distribution. Motivated by this, this sub-section is devoted to the derivation of a closed-form 
expression for the corresponding envelope pdf. It is firstly recalled that the a—p distribution is 
a generalised small-scale fading model which also accounts for the non-linearity of the wireless 
channel. It was reported in 0] and it includes as special cases the well known Weibull, Nakagami- 
m and Rayleigh distributions. 

Its envelope pdf is expressed as 


p R (r) = 




3 M rOL 


(24) 


f“PT(/i) 

The a— /i/gamma distribution is formulated by taking the conditional probability, p lt (r \ y ) and 
averaging over the statistics of the gamma distribution, namely, 


Pr(T) 


p R ( r I y)p v (y)dy 


afi^r afl 1 

Z T{n)T(b)n b 


y b-a»- i e -^ e -Hi dy 


(25) 

(26) 


B. An exact closed-form expression for the envelope pdf 

It is evident that an analytic expression for the envelope pdf of the «-/i/gamma fading model 
is subject to evaluation of the integral in (f25l) . To this end, by utilizing again (33, eq. (8.4.3.1)] 
and [f33l eq. (8.4.3.2)], it follows that, 


p R [r) = ap? 


r y1,0 (y i-) r 0 ’ 1 ( y a i 1 

u o,i (o lo) Lt 1,0 I. 


'o i/ 1+ “ /i_b r 1 - a PU fc r(/i)r(6) 


dy 


(27) 


Notably, the above integral has the same algebraic representation as (1201) . Therefore, by utilizing 
[33] eq. (2.24.1.1)], one obtains a closed-form expression for the envelope pdf of the a —///gamma 
distribution, namely, 
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To the best of the authors’ knowledge equation (l28l) is novel. 


(28) 


V. The a—n—n Extreme / GAMMA DISTRIBUTION. 

A. Model Formulation 

It is recalled that the a—n—p Extreme fading model emerges from the a-n-fi distribution and 
accounts for severe fading conditions that typically occur in enclosed environments. It constitutes 
a generalization of the k— p extreme distribution since it additionally considers the non-linearity 
effects of the wireless medium. The basic principle underlying this model is that the CLT is not 
valid since the number of multipaths is rather low. 

However, likewise the a—n — p distribution, the a—n—p Extreme fading model does not 
account for the simultaneous occurrence of shadowing. Therefore, by following the same proce¬ 
dure as above, the a—n—p/Extreme gamma distribution is formulated. Subsequently, an analytic 
expression is derived for the corresponding envelope pdf. To this end, the conditional probability 
of the ol — u —/i Extreme distribution with mode y is straightforwardly expressed as 


p R (r | y ) = 2 am 


--i 
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— 2 m—, 


yie 
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Averaging the above expression over the gamma shadowing statistics yields, 
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By expressing the Ii(.) function in terms of (fl5l) and the exponential functions according to 
eq. (8.4.3.1)] and [f33l eq. (8.4.3.2)], equation (l30l) can be equivalent re-written as 


(30) 
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Finally, with the aid of [33. eq. (2.24.1.1)], one obtains the following analytic expression for the 
envelope pdf of the a—n—y Extreme/ gamma fading model, 
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(32) 


fi. The pdf of the K—p Extreme/gamma fading model. 

The K—p Extreme /gamma composite fading model was recently reported in fT7Tl , lfl8ll . 
An alternative analytic expression for the corresponding envelope pdf can be straightforwardly 
deduced by the a—K—p Extreme /gamma distribution. This is realised for the specific case that 
a = 2. To this effect, it immediately follows that, 


PrVJ - 


E 
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2 b 2i +2 777 2 ( 1 +2/), r 2(l+0 1 T(n + l) 
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3,0 


4tt 2 


p( 1 + k)t 2 


1)1| Tj-^+M+2 \ 


-+/i+2 




(33) 



Fig. 1. Envelope pdf of the a— k— /i/gamma distribution for b = 1.1, fl = 0.9, [i = 2.1 and different values of a. 
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VI. Numerical Results and Discussions 

This Section is devoted to the demonstration of the general behaviour of the derived analytic 
expressions. Figure 1 illustrates the envelope pdf of a— k— /i/gamma distribution with respect to 
x for the case that b = 1.1, Q = 0.9, // = 2.1 and different values of a. One can clearly observe 
the flexibility offered by the non-linearity parameter. Likewise, Figure 2 considers b = 1.8, 
0 = 0.7, k = 4.0, a = 2 and different values of //. One again can observe the sensitivity of the 
fi parameter particularly in the small value regime. 

In the same context, Figure 3 demonstrates the envelope pdf of the a— /i/gamma distribution 
for the general case that b = 1.1, 0 = 0.9, // = 2.1 and different values of a. The characterization 
flexibility thanks to a is evident and this model is considered ideal for small-scale fading in 
NLOS communication scenarios. Finally, the behaviour of the envelope pdf of the a—K—[i 
extreme/ gamma fading model is illustrated in Figure 4 for Q = 0.8, b = 1.2, m = 1.1 and 
different values of a. 



Fig. 2. Envelope pdf of the re—/i/gamma distribution for a = 2, b = 1.8, U = 0.7, re = 4.0 and different values of fi 
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Fig. 3. Envelope pdf of the a—fi fading model for b = 1.1, SI = 0.9, ^ = 2.1 and different values of a 

A. Usefulness and Applicability in Wireless Communications 

It is widely known that the algebraic representation of crucial performance measures is rather 
critical in studies related to analytical performance evaluation of digital communications. This 
is obvious by the fact that when the algebraic form of a corresponding measure is convenient, 
it ultimately becomes more possible that the derived relationships can be expressed in closed- 
form. Based on this, the fact that the form of the offered analytic expressions is suitable due to 
the existence of various identities involving Meijer G-functions, constitutes the proposed fading 
models convenient to handle both analytically and numerically. Therefore, the derived expressions 
can be efficiently applied in various analytic studies relating to the performance evaluation 
of digital communications over composite multipath/shadowing fading channels including non- 
linearities as well as severe fading conditions. Indicatively, they can be meaningfully utilized 
in deriving explicit expressions for important performance measures such as, error probability, 
probability of outage, ergodic capacity, channel capacity under different transmission policies and 
higher order statistics. It is recalled here that expressions corresponding to the aforementioned 
measures can obviously be derived in both classical and emerging technologies such as single 
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channel and multichannel communications, cognitive radio and cooperative systems and free- 
space optical communications, to name a few. 



Fig. 4. Envelope pdf of the a—K—fi /gamma distribution for Q. = 0.8, b = 1.2, m = 1.1 and different values of a. 

VII. Closing Remarks 

This work was devoted to the introduction, formulation and derivation of the a—K—mu /gamma 
composite fading distribution. By using this as a basis, the a— //i/gamma and the a—K—fi 
extreme /gamma distributions where allso formulated and derived. These models assume gamma 
distributed shadowing and are particularly flexible including as special cases the Weibull/gamma, 
Nakagami-m /gamma, Rayleigh/gamma and Rice/gamma multipath/shadowing fading models. 
Novel analytic expressions were derived for the envelope probability density function which 
can be considered useful in applications related analytical performance evaluation of digital 
communications over generalized multipath/shadowing environments with non-linearities. 
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